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Abstract 

We present a possibility to switch the chirality of a spin vortex occurring in a magnetic nanodisk 
by applying a uniform in-plane field pulse, based on optimizing its strength and duration. The re- 
lated spin-dynamical process, investigated by micromagnetic simulations, consists of several stages. 
After applying the field, the original vortex is expelled from the disk, after which two C-shaped 
states oscillate between each other. The essence of the method is based on turning the field off at 
a suitably chosen moment for which the orientation of the C-state will evolve into the nucleation 
of a vortex with the desirable chirality. This idea simply uses the information about the original 
chirality present inside the nanodisk during the dynamic process before losing it in saturation, and 
can thus be regarded as analogous to the recent studies on the polarity switching. 
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Recent progress in technologies of preparing confined nanomagnetic geometries^* 2 - has 
considerably improved the parameters of spin-related devices and shifts the challenges in 
nanomagnetism towards new approaches.-^ Several interesting possibilities emerged from 
extensive research on spin vortices occurring at equilibrium of flat ferromagnetic cylinders 
and analogous elements,- 16 - whose two binary characteristics, the chirality c = ±1 (counter- 
clockwise (CCW) or clockwise (CW) vortex's flow) and the polarity p = ±1 (up or down 
orientation of the vortex's out-of-plane polarized core,- 1 ^) determine the dynamic response 
to ultrafast magnetic field pulses 9 and are important candidates for nonvolatile magnetic 
memory^ and other applications.— 1 ^ 

Switching these states traditionally involved a strong quasistatic magnetic field, out-of- 
plane for the polarity^ and in-plane for the chirality. To switch the latter, however, required 
a geometric asymmetry, because a saturated symmetric cylinder is — after turning the field 
off — ordered randomly— dkd!L On the other hand, ultrafast dynamic processes, being forced 
by field pulses^^^ 9 - or electrical currents,— i2ii22i22i2i enable to overcome those problems by 
utilizing consecutively occurring stages. Thus the polarity can be switched via the generation 
of a vortex-antivortex pair by a short in-plane field pulse of a relatively low amplitude, 
because the original vortex then annihilates with the antivortex, leaving there the new 
vortex with the opposite polarity. 

We here show, using micromagnetic simulations, that the chirality of a vortex in a sym- 
metric nanodisk can also be switched by an in-plane field pulse, because the information 
about the original state is still present inside the disk for a certain time (before arriving at 
the saturated equilibrium), so that a vortex with the opposite chirality can be nucleated 
after turning the field off at a suitably chosen moment. This method eliminates the need 
for inhomogeneity of the excitation, which has been introduced by a mask, 25 spin transfer 
torque, 26 or exchange bias.— 

We carry out the time integration of the Landau-Lifshitz-Gilbert (LLG) equation 

dm dm 

—— = - 7m x // + am x — — , (1) 

at at 

where m denotes the unit magnetization vector, 7 the gyromagnetic ratio, a the Gilbert 
damping parameter, t the time, and H c s = —(fioM s ) dE /dm the effective magnetic field, 
determined from the energy density E (with p, Q being the magnetic permeability of vac- 
uum and M s the saturation constant) and being the sum of the externally applied field 
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term H ex t, exchange term i/ CX ch = {2A/ fi( ) M s )'V 2 m (with A denoting the exchange stiff- 
ness constant), and demagnetization term The LLG equation and all the field terms 
are numerically treated via regular discretization.— Thus the exchange field is calculated as 
-Hcxch = (2A/3[ioM s d 2 ) Trfi\ with the summation over all eight neighbours in the grid.— 
The demagnetizing field is generated by surface magnetic charges present on each cell, with 
help of analytical formulae used for surface integrations as described by Hubert and Schafer.- 
The method was tested to be very close to the OOMMF public code.— All numerical exper- 
iments were carried out on a Py nanodisk of the diameter of 200 nm, thickness of 20 nm, 
with c = p = 1, using material parameters M s = 860 kAm -1 , A = 1.3 x 10 -11 Jm _1 , and 
a = 0.01. 

Trajectories plotted in Fig. QJa) demonstrate the initiation of motion of the vortex core 
in the Py nanodisk according to the amplitude of an applied in-plane field pulse. For fields 
below 50.6 mT, the vortex core starts to move along the field and then traces a CCW 
spiral around a new equilibrium point, shifted from the disk's center (45 mT chosen as an 
example). For fields between 50.6 and 57.7 mT, the spiral motion leads the vortex core out 
of the disk, so that the vortex annihilates. For higher fields a new vortex of the opposite 
polarity and an antivortex are created, the latter of which annihilates with the old vortex, 
facilitating thus the switching of the polarity. According to the classification in Fig. [T](b), 
to switch the chirality obviously requires intermediate fields (examples in Fig. [Ha) chosen 
51 and 57 mT). 

The whole process after applying a finite field pulse with the amplitude of 51 mT is 
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FIG. 1: (color online) Selected trajectories of the vortex core plotted inside the top- viewed contour 
of the Py nanodisk (a) for applied fields B x = 45, 51, and 57 mT. For the case of 45 mT (dashed 
curve) the trajectory corresponds to a spiral motion, whereas for 51 and 57 mT (solid and dotted 
curves) the trajectory leads the vortex core out of the disk. The chart (b) classifies all processes 
according to the size of the field. 




3 




FIG. 2: (color online) Time evolution of magnetization distribution after applying a pulse of 
magnetic field B x = 51 mT with two examples of temporal duration of 836 and 988 ps. The color 
scale of M y , the geometric orientation, and the field pulse are displayed in the bottom- left corner. 
The top part of the figure demonstrates the evolution during the nonzero field, applied from t = 
until a stable, saturated state at 5 ns. The blue (gray solid) frame displays the evolution after 
turning the field off at t Q g = 836 ps (shorter pulse). The red (gray dashed) frame displays the 
evolution after turning the field off at t' oS = 988 ps (longer pulse). 

demonstrated in Fig. [2J with two examples of pulse durations of 836 and 988 ps. As already 
shown in Fig. [TJ the CCW vortex follows a curve towards the annihilation point at 584 ps, 
after which a CCW C-like curved state (C-state) appears, which is later transformed into 
an opposite (CW) C-state, both of which oscillate between each other. According to the 
snapshots at 682, 743, and 836 ps (and movie l.mov^), the oscillation can be described as 
alternation between the CCW and CW C-states, the first of which turns into the second via 
domain wall motion from the north pole to the south pole of the nanodisk. The curvatures 
of the two C-states continuously decrease, both finally converging into a stable, saturated 
state (displayed for 5 ns). If we turn the field off at t s = 836 ps, where the most distinct 
CW orientation of the C-state is present, then a new vortex nucleates at around t g + 305 ps 
with the CW chirality, so that the chirality of the original vortex is switched (the solid frame 
of Fig. [2J and movie2.mov2i). On the other hand, if we turn the field off at t Q s = 988 ps, 
at which the C-state possesses the CCW orientation, the original chirality is preserved (the 
dashed frame and movie3.mov^i). In other words, the orientation of the C-state at which 
the field is turned off determines the chirality of the nucleated vortex. 

The dynamic evolution can be also described by a quantity called fluid-like vorticity, 
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FIG. 3: (color online) Analyzing the dynamic processes, (a) Time evolution of fluid-like vorticity 
<jj = [V x m] z (averaged over the nanodisk) after applying the pulse B x = 51 mT (the u curve) 
and after two examples of turning the field off at t g = 836 ps (the wa curve) and t' oS = 988 ps 
(the wb curve). The top inset (b) shows the Fourier transform of a selected range of the vorticity 
function u(t); the selection is depicted by the quotation arrows. The space distribution of magne- 
tization oscillations during the nonzero field for the peak frequency f^o are displayed in the bottom 
insets corresponding to the amplitudes (c) and the phases (d) of the oscillations. 

uj = [V x m] z (with the z axis oriented out-of-plane) , whose curves (averaged over the 
nanodisk) are plotted in Fig. [3J (a). Here the positive values represent spin distributions for 
which the CCW orientation predominates, whereas the negative values represent the CW 
orientation. Thus the curve u(t) exhibits a steep decrease of the vorticity for t < 836 ps, 
signifying the disappearance of the CCW vortex and evolution into the CW C-state, and 
then damped harmonic oscillation signifying the alternation between the CW and CCW 
C-states. To describe this oscillation in detail, we performed the Fourier transform of a 
selected part of the vorticity function u(t) [Fig. [3J (b)], showing a strong frequency peak 
at Qo = 2.8 GHz, broadened due to the damping and the space inhomogeneities of the 
oscillation. For the chosen frequency Qq, the spin-dynamical motion in the x-y plane can 
be mathematically described as 



where mo is the immobile component of magnetization, a denotes the distribution of am- 



m(x, y; t) = mo(x, y) + a(x, y) cos[fi ^ + S(x, y)] 



(2) 
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FIG. 4: (color online) Time evolution of energies, (a) The evolution of the total, exchange, de- 
magnetizing, and external-field energies during the externally applied field B x = 51 mT. The two 
quotation arrows show two examples of turning the field off. (b) The evolution of the energies after 
turning the field off at t ^ = 836 ps. (c) The analogous evolution for t' oS = 988 ps. The trajectories 
of the cores of nucleated vortices, corresponding to the two cases (b) and (c), are displayed in (d) 
and (e), respectively. 

plitudes and 5 the distribution of phases of the oscillations within the nanodisk. These a 
and S distributions are displayed in Fig. [3] (c, d), respectively. The highest amplitudes of 
the oscillations are obviously on the left and right edges of the nanodisk (on the left more 
significant), while the phases are constant on vertical lines. In other words, the motion is 
comparable to the napping of two wings, with a certain asymmetry of the amplitudes and 
phases due to the asymmetric origin of the vortex annihilation. These modes have been 
studied in detail previously.— 1 ^ 

For further illustration of the spin evolution after the two pulse durations of 836 and 
988 ps (the two bottom frames of Fig. [2]), we also show the corresponding evolutions of the 
vorticities, uj\(t) and u B (t), plotted in Fig. [3] (a), exhibiting the nucleation of the CW or 
CCW vortex, respectively. Both curves clearly confirm that the negative or positive sign of 
the vorticity of the alternating C-states uniquely determines the CW or CCW chirality of 
the final vortex, so that the entire process is controllable by choosing the appropriate pulse 
length, analogously to the current-driven domain wall motion studied by Thomas et al.— 

We have also investigated the reliability of the method with respect to the edge rough- 
ness^ by utilizing the deviation from the cylindrical symmetry due to the staircase approx- 
imation.—^ Thus, besides the excitation field oriented along the x-axis, we also performed 
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analogous simulations with fields rotated by 30° and 45° to reveal that the initial process 
of the vortex expulsion remains nearly same, followed by slightly faster C-state oscillations. 
However, the time t oS discussed above remains quite adequate for the chirality switching, 
regardless of the field orientation. 

To complete our demonstration, Fig. 0] displays the micromagnetic energies of the nan- 
odisk, corresponding to all the processes discussed above. Fig. 0] (a) shows the energies 
during the externally applied field pulse of 51 mT. After applying the field at t — 0, the 
vortex gradually changes its distribution to reduce the external field energy, while the de- 
magnetizing energy increases. The exchange energy, which is mostly coupled with the vortex 
core, significantly vanishes during the vortex annihilation. The total energy exhibits grad- 
ual decrease due to damping. Fig. [4] (b) shows the energies after turning the field off at 
t s = 836 ps, exhibiting a rapid increase of the exchange energy and a decrease of the 
demagnetizing energy because a new vortex nucleates. The two red asterisks denote two 
points of strong release of energy due to the switching of the polarity of the new vortex, 
accompanied with the creation and annihilation of vortex-antivortex pairs. Fig. H] (c) shows 
analogous processes for the case t' oS = 988 ps. Fig. 0] (d, e) displays the trajectories of the 
created vortex cores for the two cases corresponding to Fig. 0] (b, c), respectively. 

In summary, we have proposed a switching method for the vortex chirality based on op- 
timizing the strength and duration of the excitation pulse. The method does not require 
any artificial asymmetry of the nanodisk or any particular field distribution, because the 
original chirality determines the entire process. The idea is particularly appropriate for nan- 
odisks with relatively small diameters where easy vortex expulsion is possible without core 
switching. For large diameters, however, the method will also require optimizing the pulse 
shape or combining the dynamic pulse with quasistatic field. A weak point of the method 
is the fact that the information about the original polarity is lost after the annihilation of 
the original vortex, which can be prevented by memorizing the polarity before the vortex 
is excited and (eventually) by switching it after the new vortex relaxes. Another problem 
might arise from the fact that the excitation pulse must be applied with a relatively precise 
duration, which can slightly differ for different shapes of the nanodisk occurring due to fab- 
rication imperfections; however, the minor influence of edge roughness does not cause much 
trouble. 

This work is part of the research plan MSM 0021620834 financed by the Ministry of 
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